Spontaneous Spin Polarization in Quantum Wires 
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A number of recent experiments report spin polarization in quantum wires in the absence of 
magnetic fields. These observations are in apparent contradiction with the Lieb-Mattis theorem, 
which forbids spontaneous spin polarization in one dimension. We show that sufBciently strong 
interactions between electrons induce deviations from the strictly one-dimensional geometry and 
indeed give rise to a ferromagnetic ground state in a certain range of electron densities. 

PACS numbers: 73.21. Hb,73.63.Nm,75.10.Pq,75.30.Et 



Quantum wires are quasi-one-dimensional structures 
which, although conceptually simple, display extremely 
rich physics that defies conventional intuition developed 
for two- and three-dimensional conductors. The study of 
transport properties of quantum wires has offered a di- 
rect glimpse into the quantum world through the quanti- 
zation of conductance in integer multiples of Go — 2e^ /h 
Q. Recently, one of the most exotic implications of 
one-dimensionality — the existence of separate spin and 
charge excitations — has been demonstrated experimen- 
tally 0. 

In a number of recent experiments on quantum wires, 
deviations frornperfect conductance quantization have 
been observed [1 HIE S S S 0- Most commonly the 
experimental findings have been interpreted as indication 
of spontaneous spin polarization 0, 0, S IE S ■ How- 
ever, for a strictly one-dimensional system this possibil- 
ity is explicitly forbidden by a theorem due to E. Lieb 
and D. Mattis based on very general mathematical 
properties of the Schrodinger equation describing these 
interacting electronic systems. Although a number of in- 
terpretations of the conductance anomalies that do not 
rel y on the idea of spin polarization have been proposed 
[3, lul llJ) ll^ ) the experiments do raise a fundamental 
question: Can the ground state of the electron system in 
a quantum wire be ferromagnetic ? 

The only way to circumvent the Lieb-Mattis theorem 
is to recognize that realistic quantum wires are not in 
essence one-dimensional devices. Attempts in that di- 
rection have been made [l^ . requiring, however, a fully 
two-dimensional structure as a starting point. By con- 
trast, we start with the conventional model of an electron 
gas in a quantum wire and show that strong Coulomb in- 
teractions both cause deviations from one-dimensionality 
and bring about a ferromagnetic ground state. 

Typical experiments are done with quantum wires 
that are formed at the interfaces of GaAs/AlGaAs het- 
erostructures. A voltage applied to metal gates provides 
a confining potential in the directions transverse to the 
wire and, in addition, allows one to tune the electron den- 
sity in the wire. While conductance plateaus at integer 
multiples of Go are observed in the high density reg;iniej_a 
drop in conductance commonly attributed 



to spin polarization has been observed S S IE S IE 
in the region of gate voltages where the electron density 
is very low. 

As the density n of electrons is lowered. Coulomb in- 
teractions become more important, and at n ^ a^^ they 
dominate over the kinetic energy. (Here ob = h?t/me^ 
is the Bohr radius in the material, e is its dielectric con- 
stant, and m is the effective electron mass; as « lOOA m 
GaAs.) In this limit the electrons can be viewed as classi- 
cal particles. In order to minimize their mutual Coulomb 
repulsion, electrons occupy equidistant positions along 
the wire, forming a structure with short-range crystalline 
order — the so-called Wigner crystal. Upon increasing 
the density, the inter-electron distance diminishes, and 
the resulting stronger electron repulsion eventually over- 
comes the confining potential, transforming the classical 
one-dimensional Wigner crystal into a staggered or zig- 
zag chain [lE|. Typical structures for different densities 
are shown in Fig. 1. 

Quantum-mechanically, spin-spin interactions in the 
Wigner crystal arise due to exchange processes, in which 
two electrons switch positions by tunneling through the 
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FIG. 1: Wigner crystal of electrons in a quantum wire defined 
by gates (shaded) . The structure is determined by the param- 
eter V proportional to electron density (see text). As density 
grows, the one-dimensional crystal (a) gives way to a zig-zag 
chain (b-d). The arrows in (b) illustrate the nearest neighbor 
(Ji) and next-nearest neighbor (J2) exchange processes. 
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potential barrier that separates them. The barrier is 
created by the two exchanging particles as well as all 
other electrons in the wire. Originating in tunneling, the 
exchange energy associated with such processes falls off 
exponentially with the distance between the electrons. 
As a result, only the nearest-neighbor exchange is rele- 
vant in a one-dimensional crystal. The corresponding ex- 
change constant is positive, leading to an antiferromag- 
netic ground state in accordance with the Lieb-Mattis 
theorem [foj . 

A very different situation arises when one considers 
the most trivial deviation from the one-dimensional crys- 
tal, namely the zig-zag chain introduced above. For that 
structure, depending on the distance between the two 
rows which varies as a function of density, the distance 
between next-nearest neighbors may be equal to or even 
smaller than the distance between nearest neighbors, as 
illustrated in Fig. l(c,d). Accordingly, the next-nearest 
neighbor exchange constant J2 may be equal to or larger 
than the nearest neighbor exchange constant Ji. The 
corresponding spin chain is described by the Hamiltonian 



H12 - ^ (JiSjSj+i 



J2S1S 



3^3+2) 



(1) 



The competition between the two exchanges causes frus- 
tration of the antiferromagnetic spin order and eventually 
leads to a g apped dimerized ground state at J2 > 0.24Ji, 
[idL ITtI |l8i| . In addition, drawing intuition from studies 
of the two-dimensional Wigner crystal, one realizes that 
in this geometry ring-exchange processes, in which three 
or more particles exchange positions in a cyclic fashion, 
have to be considered. 

It has been established that, due to symmetry proper- 
ties of the ground state wave functions, ring exchanges 
of an even number of fcrmions favor antiferromagnetism, 
while those of an odd number of fcrmions favor ferromag- 
netism In a zig-zag chain, the Hamiltonian reads 

^ = 2 X! (^l^J J+l + -^S-Pj j-h2 - J3{Pj3+i3+-^ + Pj+'ij+ij) 
3 

-I- J4(Pj j+l j+3 + P3+2j+3j+lj) (2) 

where the exchange constants are defined such that all 
Ji > and only the dominant Z-particle exchanges are 
shown. Here, Pji...ji denotes the cyclic permutation op- 
erator of I spins. A more familiar form of the Hamil- 
tonian in terms of spin operators is obtained using that 

Pij ~ 2 ~^ 2SiSj and Pji...ji — Pjij2P3233 ■ ■ ■ P31-131 E3- 
In particular, the two-spin exchanges reduce to Eq. 

The simplest ring exchange involves three particles and 
is therefore ferromagnetic. Extensive studies of the two- 
dimensional Wigner crystal have shown that, at low den- 
sities (or strong interactions), the three-particle ring ex- 
change dominates over the two-particle exchange. As 



a result, the two-dimensional Wigner crystal becomes 
ferromagnetic at sufficiently strong interactions [20I l2l| . 
Since the electrons in a two-dimensional Wigner crystal 
form a triangular lattice, by analogy, one should expect 
a similar effect in the zig-zag chain at densities where 
the electrons form approximately equilateral triangles, 
Fig. 1(c). In order to verify this scenario, we have to 
identify the electron configuration that is stable at a 
given density and subsequently find the corresponding 
exchange energies. 

Specifically, we consider a quantum wire with a 
parabolic confining potential Vconf (y) = where 
n is the frequency of harmonic oscillations in the po- 
tential Vconf(2/)- At low electron density n in the wire, 
a one-dimensional Wigner crystal is formed. Fig. 1(a). 
As the density grows, however, the Coulomb interaction 
energy becomes comparable to the confining potential, 
leading to the formation of a zig-zag chain, as depicted 
in Fig. l(b-d). This transition happens when distances 
between electrons are of the order of the characteristic 

1 /3 

length scale ro — (2e^/eTOil^) , such that K:onf(?'o) = 
Mnt(''o)j where Vint{r) = e^/er is the Coulomb interac- 
tion energy. It is convenient for the following discussion 
to introduce a dimensionless density v — nro. Minimiza- 
tion of the energy with respect to the electron configura- 
tion reveals that a one-dimensional crystal is stable 
for densities v < 0.78, whereas a zig-zag chain forms 
at intermediate densities 0.78 < v < 1.75. (At higher 
densities, the zig-zag chain gives way to structures with 
larger numbers of rows 0|.) The distance between rows 
grows with density, and the equilateral configuration is 
achieved at v k, 1.46, well within the region where the 
zig-zag chain is stable. Therefore, there are strong indica- 
tions that the ferromagnetic state may be realized. More 
specifically, upon increasing the density one would ex- 
pect the system to undergo two consecutive phase transi- 
tions: first from an antiferromagnetic to a ferromagnetic, 
and then to a dimer phase. However, the latter scenario 
cannot be established conclusively based solely on the 
two-dimensional Wigner crystal physics. The main dif- 
ferences are (i) the presence of a confining potential as 
opposed to the fiat background in the two-dimensional 
case, and even more importantly, (ii) the change of the 
electron configuration with density. Fig. 1, as opposed to 
the ideal triangular lattice in two dimensions. Below, we 
study numerically the exchange energies for the specific 
configurations of the zig-zag Wigner crystal in a parabolic 
confining potential. 

The strength of the interactions is characterized by the 
parameter 



_ ro_ _ / me^ 1 



2/3 



(3) 



For rn ^ 1, the physics of the system is dominated 
by strong interactions, and a semiclassical description is 
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applicable. In order to calculate the various exchange 
constants, we use the standard instanton method, also 
employed in the study of the two-dimensional Wigner 
crystal li^,!!^- Within this approach, the exchange con- 
stants are given by J/ = J* exp (—Si/h), where Si is the 
value of the Euclidean (imaginary time) action, evaluated 
along the classical exchange path. By measuring length 
and time in units of rp and T = \/2/^, respectively, the 
action 5'[{rj(T}] is rewritten in the form S — hrj^Jrn, 
where the functional 



(4) 
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is dimensionless. 

Thus, we find the exchange constants in the form Ji — 
J* exp {—rji^/rn), where the dimensionless coefficients rii 
depend only on the electron configuration (cf. Fig. 1) or, 
equivalcntly, on density v. The instanton trajectories, 
and subsequently the exponents r]i, are calculated for 
each type of exchange by solving the equations of motion 
obtained from the dimensionless action Q numerically. 

To first approximation, we neglect the motion of all 
"spectators" — the electrons in the crystal to the left and 
to the right of the exchanging particles. Figure 2 shows 
the calculated exponents for various exchanges as a func- 
tion of dimensionless density z/. At strong interactions 
(rn ^ 1), the exchange with the smallest value of ry; is 
clearly dominant, and the prefactor J;* is of secondary 
importance to our argument. The numerical calcula- 
tion confirms our original expectation: the dominant ex- 
change constant changes from nearest neighbor exchange 
Ji to three-particle ring exchange J3 to next-nearest 
neighbor exchange J2. More complicated ring exchanges 
have also been computed. Figure 2 displays the ones with 
the smallest exponents, namely the four-particle ring ex- 
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FIG. 2; The exponents rji as functions of the dimensionless 
density computed with frozen spectators. The insets illus- 
trate the four most important exchange processes. 



FIG. 3: The phase diagram including nearest neighbor, next- 
nearest neighbor, and three-particle ring exchanges. The ef- 
fective couplings Ji and J2 are defined in the text. The solid 
line shows schematically the traversal of the various phases 
with increasing dimensionless density as dictated by the 
calculated exchange energies. 



change as well as five-, six-, and seven-particle ring ex- 
changes (dashed lines). 

If one includes only the dominant exchanges Ji, J2, 
and J3, the Hamiltonian of the corresponding spin chain 
takes a simple form. Nearest and next-nearest neighbor 
exchanges are described by Eq. Furthermore, the 

three-particle ring exchange does not introduce a new 
type of coupling, but modifies the two-particle exchange 
constants [l9j. For a zig-zag crystal we find 



H3 - -Js^ (2SjSj+i + SjS 



i+2j 



(5) 



Thus the total Hamiltonian still has the form (Q, but 
with the effective two-particle exchange constants Ji = 
Ji — 2J3 and J2 = J2 — Ja- Therefore, the regions of 
negative (i.e. ferromagnetic) nearest and/or next-nearest 
neighbor coupling become accessible. The phase dia- 
gram of the Heisenberg spin chain with both posi- 
tive and negative couplings is well studied 0, 0, 



l23l |2J, l25l |26j . In addition to the antiferromagnetic and 
dimer phases discussed earlier, a ferromagnetic phase ex- 
ists for Ji < max{0, — 4J2} The phase diagram 
in terms of the effective exchange constants Ji and J2 
is shown in Fig. 3. The solid line represents schemati- 
cally the path followed in phase space, according to our 
numerical calculation of the exchange constants, as the 
density ly increases. At low densities, the system is close 
to one-dimensional and is, therefore, antiferromagnetic. 
In the range of densities corresponding to an "approxi- 
mately equilateral" configuration, the three-particle ring 
exchange is strong, leading to a ferromagnetic ground 
state. Finally, at even higher densities, frustration caused 
by the next-nearest neighbor coupling J2 drives the sys- 
tem into a dimerized phase. (Note that there is some con- 
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FIG. 4: The exponents rji, r]2, rjz, and 774 as functions of the 
dimensionless density v. The computation includes 12 moving 
spectator particles on either side of the exchanging particles. 
Corrections to rji from the remaining spectators do not exceed 
0.1%. 

troversy concerning the physics of the parameter regime 
— 4 J2 < Ji < 0, where the existence of a spectral gap as- 
sociated with dimerization has not yet been established 
conclusively 

It turns out that the above picture, based on the cal- 
culation of the exponents to first approximation, is in- 
complete: because only the exchanging particles were al- 
lowed to move while all spectators were frozen in place, 
the values of rji were overestimated. Surprisingly, allow- 
ing spectators to move results not only in quantitative 
but in qualitative changes as seen in Fig. 4. At large 
densities, the four-particle ring exchange J4 dominates 
over J2. Contrary to J3, the four-particle ring exchange 
not only modifies the nearest and next-nearest neighbor 
exchange constants — in addition, it introduces more com- 
plicated spin interactions |l9l |. For the zig-zag chain, we 
find 

3 ^ ^ 

Hi = SjSj-f; + 2[(SjSj+i)(Sj^2Sj+3) 

+ (SjS^2)(S^iS^3) - (SjSj+3)(S^iSj+2)])- (6) 

Not much is known about the physics of zig-zag spin 
chains with interactions of this type. Preliminary nu- 
merical studies indicate that the ground state has zero 
magnetization ^7]. Further work is required to identify 
the possibly novel spin structures. We would also like 
to point out that a confining potential of different shape 
might alter the outcome of the competition between the 
very close values of 774 and 772 at high densities. 

In experiments with quantum wires, the interaction 
strength is not a tunable parameter: it is determined 
by the electron charge e and the dielectric constant e in 
the semiconductor host. However, the parameter can 
still be tuned by adjusting the confining potential. As 



rn oc fi"^/'^, making the confining potential more shallow 
effectively increases interaction effects. Quantum wires 
in semiconductor heterostructures are fabricated using 
either cleaved-edge-overgrowth or split-gate techniques. 
In cleaved-edge-overgrowth wires |3, we estimate that 
rji is at most of order unity due to the steep confining 
potential. A more shallow confining potential is achieved 
in split-gate wires 0, S IE IE ■ Using the device specifi- 
cations of Ref . 13 , one obtains values of rn in the range 
ro « 3 — 6. It is not clear whether these values are large 
enough to result in spontaneous spin polarization. The 
ideal devices for observation of ferromagnetism would be 
ultra-clean wires with widely separated gates to provide 
the most shallow confining potential possible. 

In conclusion, interactions lead to deviations from one- 
dimensionality in realistic quantum wires and, as a con- 
sequence, the Lieb-Mattis theorem no longer applies. We 
have shown that strong enough interactions induce a fer- 
romagnetic ground state in a certain range of electron 
densities, where the electrons form a zig-zag Wigner crys- 
tal. 
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